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The stress distribution in an incomplete tore loaded as shown in 
Fig. 1 is of particular interest since it very closely approximates that 
in heavy close-coiled helical springs under axial tension or compression, 
Necessarily the spring helix angle must be small, which is the case in a 
close-coiled spring. By a heevy spring is meant one whose ratio of mean 
diameter to cross-sectional diameter is such a value that the curvature 
of the section must be considered. 

It should be noted that the 
stress distribution arising from 
the loading in Fig 1. is not pure 
torsion in the usual sense, but 
is a combination of torsion and 
direct shear. The problem there- 
fore resolves itself into one of 
finding a single stress function which defines the true stress distribution 
in a cross-section of the circular ring sector. 

Several solutions to the problem are in the literature, all of which 
by various means solve the differential equation arising from the conditions 
of compatibility. The first, by Hlchell (1) in 1699, used polynomial 
stress functions and obtained solutions for approximately circular cross- 
sections. Gbhner (2) used successive approximations to approach an exact 
solution. Shepherd (3) used a method similar to both Gdhner and Mlchell 
by finding a sequence of functions for approximately circular cross-sections 
and combining them linearly in such a manner that the sum was a solution. 




Fig. 1. 
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Wahl (AJ obtained a solution using curved bar theory and assuming a dis- 
placement of the center of rotation . Southwell (5) presented a formal sol- 
ution for an arbitrary cross-section with a view towards a "relaxation M 
approach. Prieberger (6) has presented an exact solution for a circular 
cross-section by finding a stress function analogous to the ordinary torsion 
function and solving the problem in toroidal harmonics. 

In this paper an approximate solution is obi&inod using the principle 
of least work. A stress function is found satisfying the equations of equil- 
ibrium and the boundary conditions and whose corresponding stresses make the 
strain energy a minimum. The solution of the differential equation of compat- 
ibility has therefore been replaced by the problem of minimising the strain 
energy. In the energy method , the condition of minimum strain energy is 
equivalent to satisfying compatibility not in a point by point sense, but 
"on the average" throughout the body. 

The purpose of this investigation has bsen to answer two questions in 
the author's mind. Namely, in view of the fact that nowhere wac the author 
able to find the energy method used in the literature: 

(1) Can the problem be solved by this method, and how do the results 
compare with those of other solutions? 

(2) Coes the problem particularly lend itself to solution by 
energy methods? 

It was f out id that the problem is not adaptable to an exact solution by 
energy methods, but by making some approximations, excellent results are 
obtained that agree very closely with Frieberger's exact solution. 
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FORMULATION OF THE PROBLEM 



We will consider ft sector of a circular ring with mean radius of curva- 
ture E and cross-sectional radius a. A load P is applied to one terminal 
cross-section as shown in Fig. 2, the other remaining fixed. Cylindrical 
coordinates are used, where the * axis coincides with the toroidal axis. 




Fig. 2. 



6 increases positively as shewn in the figure and r increases outward 
from the toroidal axis. Later in the solution the coordinates will be trans- 
formed, but for the present purpose of establishing a stress function satis- 
fying ths equations of equilibrium, cylindrical coordinates are most convenient 
Assuming zero body forces, from THEORY OF ELASTICITY, Timoshenko & Qoodier 
Equations (170) the differential equations of equilibrium are 
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Using the same assumptions mads by Gdhner in this ease, namely that the 
enly non-vanishing stresses are 1 T 2fi and that the stress distribution in 
any cross-section is independent of © these reduce to 

3Z r ° 

This siay also be written 

^ (^T xe jj * 0 

A stress function satisfying the above is 

GR 1 - || ■ <r-T^ g GR" |£ > - ^ T is 

Where £ is a constant (actually the modulus cf rigidity). 

Therefore the stresses may be expressed as 
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At this point it is convenient to transform the cylindrical coordinates 
££§ into toroidal coordinates ^ ^4^ & (refer to Fig. 3). 
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In a plane croes-scction determined by 0 a constant 
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Using liquations (2), (3) and (A) the following result is obtained. 
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The latter expressions relate the stress function and the stresses in 
the new system of coordinates « 

It follows that since the shear stress / Tj 9 0 is normal to the boundary, 
it must vanish everywhere on the boundary. This is true because the surface 

of the body is free from any external forces. Using this condition with 

94 * 

Equation (5), it is apparent that ^ = O and <P must be constant on the 
boundary. 

The circular ring sector we are considering is a sirigly connected body, 
hence the constant may be chosen arbitrarily. Therefore the boundary condition 
is taken as ^ r O everywhere on the boundary. 

The only action ©n a cross-section Is a force P directed along the 
toroidal axis. This may be resolved into a force and a couple as shown in 
Fig. 4. 




It is now seen that the two conditions of static equilibrium to bo 
satisfied are that the resultant stress on a cross-section produce & force P 
directed along the z axis end a moment about the center PH. These requirements 
may be written * X v 

pe s ir\ f Tlj, $ cos p ^ 'P 

P9 p* ± 
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The strain energy P«* unit angle 8 is 



a aw 

(7 ) — u- 4 // ♦ T « X«- P co ^) p ■ <f ** 

<5 C 

The method of solution will now be to take the stress function in the fora 

4>* X , wher® 4 L tiro suitably selected functions of P and 

ceo ' 

, each of which satisfies the boundary condition 4^-0 whan p = cx> . 
The coefficients L are constants which are evaluated from the minimum 
condition of strain energy. 
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FIRST APFaCXXKATION 



For a first approximation we shall take a function 4* * satisfying 
the boundary condition that it vanish everywhere on the boundary, in the 
form 4* c C cos ) . The reasons for this 

particular choice are discussed in Appendix A. Taking the partial deriva- 
tives of 4* with respect to the two variables p and 4* 

— !*■ - 



a4> 



s 2 P* c ♦ ** (^P -Q-^ COS ^ 
r R 



Substituting in Equations (5}» the following expressions are 
obtained for and T<p© • 



The appearance of the term (p_p CO s^Y' 1® the strees equations 
makes the integration required in (6) and (7) very complicated and the 
results largely unmanageable in the evaluation of the unknown coefficients 
in 4 . (See Appendix 3). Thie is particularly true when additional terms 
are used in P for a higher order of approximation, az;d in the evaluation 
of the strain energy where the stresses appear as squared terss 0 
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is always less than unity, we may write 
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Utilising this expansion, the exact stress expressions (5) may be 



approximated as follows 
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This particular fora of approximation accomplishes the desired result 
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of limiting the highest power to which the ratios and appear in the 



stress equations. 
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The partial derivatives are 
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Substituting, we arrive at the following approximate expressions for 

t w . 

Substituting these values of Tj# and T^in the first of Equations (6), 
and integrating we obtain 
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The same result is obtained from the second condition of Equations (6). 
It follows that °<4 is fixed by the requirements of static equilibrium 
and e*’ , may now be determined by the condition of minimum strain energy 
that — = O . 

Prom Equation (7) 
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Substituting the stresses from Equation (9) and integrating 
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Using these results in Equations (9) we arrive at the expressions 
for the first approximation of the stress distribution in a cross-section 
of the incomplete tore 




At the point of maximum stress where p - cx and ^ - O the above 
reduce to r 



01 ). 



i 



Tpe * ° 





It is interesting to note at this point that for this particular solution, 
one of the unknown coefficients in is determined directly from the require- 
ments of static equilibrium, and the other directly from the minimum strain 
energy condition without constraint arising from static equilibrium. 
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S&COKD APPROXIMATION 



A closer approximation to the true stress conditions will result 
if higher order terms of a suitable nature are used in the stress function* 

We shall now take £ as 

4 >« (p 1 - 0 ^ 6*0 «■ < ^cos4' v + li £ L Sm x + + A 

\ R r-*- r' R* J 

Reasons for this particular choice of functions are discussed in Appendix A. 

Again employing the binomial expansion of , — ■ — n v we write approximate 

(R-pcosvj 

expressions for 'T“p® and T • 

V* - f [(' ♦ *■*“* * -. » 0 * *?£<>.* * 4*^] 

T W ,a [0 ^w,.+ ^4 v 4 ( 

Where 

4.= Cf.J) <v 4,- 

4= 4>j= 

This is an extension of the device used before to limit the highest 
power to which the ratios and appear in each term of the stress 
equations. Since ^ and -j^ occur in a like manner in 4^ 1 4 > 3 and 4 > 4 , 
these latter terms are grouped together and treated in similar fashion when 
introduced into the approximate expressions for the stresses. 

Taking the partial derivatives, substituting and rearranging the terms 
for convenient integration, the following approximate expressions for I'p© 
and T^ e are obtained. 
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From the first of the static equilibrium conditions in (6) (that the 

resultant stress oust produce a force P in the Z direction) it is again 

PR 

found that cf * — A 

* &TTO. 

The second static equilibrium condition (that the resultant stress 
must produce a moment about the oenter equal to PR) gives the following 
result. 



PR 

GTTOl " 




However, since ^ . PR 



• (tfrv, * a t>. * • (W&) -« 



(13)' 



Since <> <*3 and ^4 will ultimately all contain the factor 

PR 

some simplification of the algebra will be afforded if we make the 

G» lT<X 

following substitutions 

(V ^ wh,r * ' > 2 j 3 

Finally the constraining function derived from the conditions of static 
equilibrium to be used in minimising the strain energy is 

(x !»)£,+ (fc I*) ^ (z V) P 3 + (?>) (^4 + fx' 



L« O 



The work involved in obtaining the partial derivatives of the strain 
energy with respect to the unknown coefficients oi i ^ } ^ and ^ ^ will 

be simplified by differentiating under the integral sign. 

Therefore ^ ^ 

sL* i //( r r^: + 

« « 7 

Sir 

— is not required since ®i 0 has already been evaluated from the 
9^0 
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conditions of static equilibrium. Since 
it is convenient here to take '‘LIT 



- (Constant) 

3(3* 

After substituting for t ; 



and 






a(3„ " V« 

from Equations (12) and performing the required integration, 



the following expressions are obtained. 
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To minimize the strain energy and evaluate the unknown coefficients 
(3i » » (3^ and (^4 * met kod Lagrangian multipliers will be 

used with the constraining function (13) established by the requirements of 

w 

static equilibrium. The constant i*~ appearing in the partial derivatives 

GR 

of the strain energy will be incorporated in the multiplier. Letting X be 
a Lagrangian multiplier, and J(j3, > {3 t ^ |3 3 , ^4^ s O c © nst raining function, 
we may write 

3U 



3^ + X 3|L 0 



where 0 ® I , z 5 3 , 



Using (13) and (14) in the above, and the fact that 
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w* arrive at the following set ©f equations, the solution of which will 
evaluate @i s Pi j j (J 4 and determine the stress function. 

‘ ( ! • ,ll‘ > i P , N a t P, ♦ ^ P4 - £ |l 
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Using these results in Equations (12), we may now write expressions 
representing a second approximation of the stress distribution in a cross- 
section of the in complete tore. 37 o? 



ircC 



t^f 



U7)- 



At the point of maximum stress, idlers p * <x and 4’ = o , the above 
reduce to 

> f9 * O 



[ r *»] r» lt . ' ♦C 4 *$S)r. + Ct _ * S )r'-] 



As is apparent from the foregoing development, further approximations 
utiliting additional terms in the streas function will result in extremely 
long and tedious calculations. This in itself is a limitation of this method. 
Therefore at this point, assuming the solution to be a rapidly converging one, 
we will stop and introduce actual values cf the ratio of cross-sectional 
radius to the mean radius cf curvature of the tore in order to compare results 
with other solutions. 
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ftSSULTS 



The distribution of shearing stress on a horizontal diameter is shown in 
Fig. 5 md the circumferential stress distribution in Fig. 6. In both cases a 
ratio cf li/a equal to 4 is used since this realises the worst condition (i.e. 
greatest curvature for a given cross- section) of any practical significance. 

The quantity S appearing as the ordinate in both curves is a dimensionless 
quantity and is equal to f since Tp@ vanishes on both a horizontal 

diameter and the periphery. A stress distribution representing pure torsion 
of a straight circular shaft is shown by a dotted line in both figures. It is 
seen that the maximum stress actually existing in the tore is considerably 
greater than that derived from ordinary torsion theory. Both curves are in 
good agreement with similar ones derived from the exact solution by Frieberger. 
Points from Frieberger' s curves appear as the small circles in Fig. 5 and 6. 

In comparing the results cf this solution with others, namely Gflhner, Wahl, 
and Frieberger, the point of maximum stress will be used as a reference with 
different values of tt/a. Table 1 gives values cf K in the expression 
Of*©] tvs**.' jq} M for the several solutions. 



Table 1. 
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Exact 


This Sc 


lution 


Other Approx. Solutions 


a 


Frieberger 


1st Approx. 


2nd Approx. 


Gdhner 


Wahl 


4 


1.376 


1.313 


1.371 


1.372 


1.400 


5 


1.293 


1.250 


1.287 


1.295 


1.310 


6 


1.237 


1.209 


1.234 


1.239 


1.252 


8 


1.171 


1.156 


1.171 


1.172 


1.184 


10 


1.136 


1.125 


1.134 


1.135 


1.145 
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Table 1. indicates that the energy method applied to this problem 
produces results which compare favorably with other solutions. It also 
appears that the solution converges rapidly, since only five terms were 
used in the stress function. 
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Fig. 5. 

Distibution of shearing stress on a horizontal 
diameter for a/R=l/4. (‘f’Oi'ir ). S=(T^)(rra )/P. 
Frieberger’s points are indicated by small circles. 
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Fig. 6. 

Circumferential stress distribution for 
a/R=l/4. S=(Tn>e) (tt® ^ )/p . Frieberger's points 
are indicated by small circles. 
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cckousioss 



In discussing any conclusions from this investigation, it would be 
appropriate to recall the two questions that prompted it. 

(1) Can the problem be solved by this method, and how do the results 
compare with these of other solutions? 

(2) Does the problem particularly lend itself to solution by energy 
methods? 

First, the method will work and acceptable results are obtained with a 
relatively few terms in the stress function. This is in itself worthy of note, 
since it allows a very complex problem to be attacked by the more elementary 
methods of mathematics. 

However, in reference tc the second question, there are limitations both 
inherent in the energy method and peculiar to this particular application, that 
strongly indicate the problem ia not especially adapted to a solution by energy 
methods. 

The energy method, except in unusual circumstances, does not provide an 
exact solution. Consequently, in the absence of an exact solution, there is no 
real basis for Judging the results. The fact also that the energy method requires 
minimising an integral, which is done only with extreme difficulty with any 
number of terms in the stress function, is a limitation to its adaptability. 

In conclusion then, it may be said that this solution has the value of 
arriving at very good results using a relatively uncomplicated stress function 
of only five terms. 
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AFrKKDIX A 



According to the principle of least work which is used in this solution, 
an exact stress function would require selecting from all functions that satisfy 
the boundary condition those which minimize the strain energy. 

Since in general this procedure is too difficult, a limited number N of 
suitable functions was selected to determine an approximate stress function. 

In choosing functions of P and in 4 * ^LtL , the 

x L '° 

first consideration was ths boundary condition = o when p » <x • 

This condition was satisfied by taking each 4^ to contain the factor ( p** 1 cl) 

Then 4* 0^“°^) 2- (f; V) 

' L*0 

With rectangular coordinates in mind, where c p Cos 4-’ and 

^ f the next logical step was to express ^ > *] 1 in a power 

series. The first six terms of such a series were considered, namely those 
involving 

1 > f > <1 > f > \ 



Since a horizontal diameter ( f| s © ) on a plane cross-section (8 is a 
constant) is an axis of symmetry for the 4 surface, 4 must be even in and 
not contain terms involving odd powers of t| • In general ip. will appear to all 
powers since ( ^ « o ) is not axis of symmetry. Therefore the remaining terms 
expressed as functions of p and 4 are 

l j pass 4 , J p <L ^/n' t 4 

The term in 4*4 while not consistent with this line of reasoning, 

appeared as a result of the binomial expansion used in approximating the stresses. 
It was extracted from Odhners solution where a like approximation was used. 
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Consequently the stress function was taken in the form 

Here p was replaced by -£ so that would in all cases be 

pD 

the product of a dimensionless number and the factor * 

In the first approximation the first two terns were used and in the second 
all five were introduced in 4 - 



23 



70 IP, 



srecrt i nl r»* .d sew nttlfonul eesii • d «;«oO 



•d s»a?o 



bnoaot ®<li ni 



il. oi birw Jsdd oa J- -^d bo© I,-®n q 9*s«t 

£9 * 

’•doel bn* ^duara -ee'.-<X inert a la daub© 1 * efid 
twjs bee*, •'icw wned ©wi Jeail »rid noideisixoi *qe du*ill i»ld nl 

. 4 s beoobotdni ei.-'w svxt lie 



e: 



APPENDIX B 



The appearance of the term ^ the exac t equations relating 

the stresses and the stress function gives rise to the occurrence of Integrals 

£_Zir , Air 

or tho typo f [ £ jej* and Jf in 

JJ 0 Cf?-pco3<P) v 6 { (R-pcoSlb)l 

evaluating the strain energy and in consideration of the conditions of static 



equilibrium. 

Taking the simplest fora of the first type, where ri^rl 



we have 
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Integrating first with respect to and setting R‘ c 
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(c.+ b cos 



. b- l <- |(ctt,coi»p ) ^ £ , 

(c + b COS b ^ 



Multiplying by 'V and integrating 

f 4 ? ^ = p , , ifi 

J <*-H> Cibcos9 bjcf 

f dt , . _k_ / sifl jfc. V _c_ 

J C + tCoS^ C l -b x \^c Pbcos^j T C^b^ J C-fr 
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SciiioJ 
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0»riT 



eoDci + j) 



. 2i> 
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Introducing the limits ° and Z.1T , this reduces to 

2ttc c= R 

(cr-tf-ft where W-p 




The other more complicated forms where r> f o and ^ ^ ° are integrable 
in finite terms by similar reduction methods, but it is apparent that the work 
becomes excessively involved. Also the results in the form just developed are 
not readily usable in evaluating the unknown coefficients in the stress function. 

In view of the foregoing, despite the fact that it was not actually necessary, 
it was expedient to approximate the stresses in such a man: er that the integration 
was simplified and the results put in a usable form. 

This device of approximating the stress equations compromised the requirement 
that the stresses satisfy the equations of equilibrium. However, it appears that, 
since the stresses do satisfy the conditions of minimum strain energy and static 
equilibria, and give satisfactory results, the compromise my be tolerated. 
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